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Recently, we showed that eigenphase shifts and eigentime delays near a resonance for a system of 
one discrete state and two continua are functionals of the Beutler-Fano formula using appropriate 
dimensionless energy units and line profile indices and identified parameters responsible for the 
avoided crossing of eigenphase shifts and eigentime delays and also parameters responsible for the 
eigentime delays due to a change in frame transformation. 

In this paper, the geometrical realization of the Beutler-Fano formulas is considered in the real 
three-dimensional Liouville space spanned by the Pauli matrices a^, cry, which are orthogonal 
in the sense that Tr(cri cTj) = 2Sij, where dynamic operators are vectors. Vectors corresponding to 
the background 5''' matrix, S matrix, and the time delay matrix Q form a spherical triangle whose 
vertex and edge angles are parameters pertaining to the frame transformations among eigenchannels 
of S, Q and eigenphase shifts of and S and the phase shift due to a resonance scattering. 
The cotangent laws of the spherical triangle yield Beutler-Fano resonance formulas appearing in 
eigenphase shifts and time delays. Duality holding for the spherical triangle explains the symmetry 
observed in the relations among parameters and provides a systematic way of defining conjugate 
dynamic parameters. The spherical triangle also shows the rule of combining the channel-channel 
couplings in the background scattering with the resonant interaction to give the avoided crossing 
interactions in the curves of eigenphase shifts as functions of energy. 

The theory developed in the previous and present papers is applied to the vibrational predissoci- 
ation of triatomic van der Waals molecules. 
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I. INTRODUCTION 

Eigenphase shifts Si {i = 1,2,. . . ,rt) of the S matrix, defined as 

S = Ue^'^U (1) 

have been utilized as a tool for analyzing the resonances . Eigenphase shifts and the corresponding eigenchannels 
are also extensively used in various forms in multichannel quantum-defect theory (MQDT) which is one of the most 
powerful theories of resonance In MQDT, the first derivatives of eigenphase shifts as functions of energy are used 
for various purposes [Q. The Lu-Fano plot is essentially a plot of the energy derivative of an eigenphase shift 
The first derivative of an eigenphase shift as a function of energy, called "a partial delay time" , has also studied as a 
relevant quantity to time delays in various fields . 

In spite of their wide use, studies of the behaviors of eigenphase shifts and their energy derivatives in the neighbor- 
hood of a resonance have not been done extensively, comparing to the studies of photo-fragmentation cross sections 
and S matrix itself. Eigenphase shifts in the multichannel system are known to show complicated behaviors near 
a resonance due to the avoided crossing between curves of eigenphase shifts along the energy ||,^ . In the previous 
work, detailed studies of the behaviors of eigenphase shifts and times delayed by collision were done for the system of 
one discrete state and two continua |^]. Comparing to the system of one discrete state and one continuum, the newly 
added one more continuum brings about the new phenomena of the avoided crossing between curves of eigenphase 
shifts and eigentimes delayed and of times delayed due to the change in frame transformation along the energy besides 
the resonance behavior. Thus the system of one discrete state and two continua provides a prototypical system for the 
study of two effects on the resonance phenomena, avoid crossings between eigenphase shifts and eigentimes delayed 
as one effect and the change in frame transformation as another one. 

Previous work showed that eigenphase shifts and eigentimes delayed due to the avoided crossing interaction and 
eigentimes delayed due to the change in frame transformation are functionals of the Beutler-Fano formula using the 
appropriate dimensionless energy units and line profile indices [Q. Parameters representing the avoided crossing of 
eigenphase shifts and eigentime delays were identified. Parameters representing the eigentime delays due to a change 
in frame transformation were shown to be described by the same parameters. With the help of new parameters, the 



1 



rotation angle 9 and rotation axis n for the S matrix exp[i(a + 9a ■ n)] were identified. The time delay matrix Q was 
shown to be given as Q = ^T^il + Pa • cr + P f ■ a) where the first term is the time delay due to the resonance, the 
second term is the one due to the avoided crossing interaction and the last term is the one due to the change in frame 
transformation. 

Though previous work found that behaviors of eigenphase shifts and eigentime delays as functions of energy follow 
Beutler-Fano's formula, it could not explain why they follow Beutler-Fano's formula. Since the system considered in 
the previous work is essentially two channel system (with one more discrete state), an analogy with the spin system 
was made and utilized but not fully exploited. One of the main purpose of the present paper is to exploit the analogy 
further. Especially, the homomorphism of the spin model with the three-dimensional rotation will be fully exploited to 
construct geometric structures made up of dynamical parameters for eigenphase shifts and eigentime delays and thus 
to derive the Beutler-Fano's formula geometrically. The geometrical realization clarifies the ambiguities in relations 
and unexplained meanings of dynamic parameters at the previous work since the geometric constructions appeal to 
our intuition and are thus easy to understand and provides means of viewing complicated relations in a simple way. 
This clarification of complicated dynamic relations through the geometrical realization is another main goal of this 
paper. 

Section II summarizes the previous results. Section III modifies the previous results a little suitable for the geomet- 
rical realization. Section IV gives a geometrical realization of the previous results. Section V connects the geometrical 
relations with the dynamical ones. Section VI applies the theory developed in Ref. ^ and the present paper to the 
vibrational predissociation of triatomic van der Waals molecules. Section VII summarizes and discusses the results. 



II. SUMMARY OF THE PREVIOUS RESULT 

Ref. examined eigenphase shifts and eigentime delays in the neighborhood of an isolated resonance for the system 
of one discrete state and two continua as a prototypical system for the study of the combined effects of the resonance 
and the indirect continuum-continuum interaction via a discrete state. The S matrix for an isolated resonance system 
is well known and, if the background S*" matrix is described by its eigenphase shifts as 5*"= C/^e"^** [/", takes the 
form @1 



5,. = J: (/>^«r + ) e"'^U^,. (2) 

where F;, F,i?o are the partial decay width of a resonance state into the Ith background eigenchannel [^, the total 
decay width J^i T;, and the resonance energy, respectively. Eq. is for the incoming wave boundary condition. The 
formula for the outgoing wave boundary condition differs from Eq. (|^) in that i is replaced by —i. 

By diagonalizing Eq. (^, eigenphase shifts S of the S matrix (—Ue^'^'^^U) for the system of one discrete state and 
two continua were obtained as 

2S± (E) = J2^" + (E) ± 5a [E) , (3) 

■i 

where Sr{E) is the well-known resonance phase shift due to the modification of the scattering wave by the quasi-bound 
state and given by — arctan(l/er) and 5a{E) is the one due to the modification of the scattering wave by the other 
wave through the indirect interaction via the quasi-bound state and was found to be given as a functional of the 
Beutler-Fano formula ||l^ 

cot^.(i?) = -cotA;,^^!^-^, (4) 



in the dimensionless energy scale defined by 



where F^ = 2^/r^/\ sin A^al and Ea =Eo + ^ cot (AF = Fi - F2, A?2 = <5i - 5^). Its form as a functional of 
the Beutler-Fano formula can be shown more explicitly by using the Beutler-Fano function 



/BF(6,g)^^^-4^, (6) 
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as 



cotSaiE) = { cot A;,y7ppj when < 

- cot A'{2 V /bF (ea , ga ) wlien €a>qa- 



The line profile index qa of the curve of Sa{E) is given by 

Ar 

2%/r\r^osA; 



9a = ■ (8) 



When Ea = qa, Sa = 7r/2 and the difference in abscissas of two eigenphase shift curves is tt/2 which is the largest 
separation of two curves when 6a is defined up to tt. Therefore, the fine profile index qa also stands for the energy of 
maximal avoidance of eigenphase shifts. En. (||) shows that the eigenphase sum Ss is given by (5e — J2i + Sr — Sj^ 
+ Sr, in conformity with Hazi's formula |ll[| . 

Let us define 5 by = U^SU^ . Let S be diagonalized by the V matrix composed of eigenvectors corresponding to 
(5+ and ^_ as 1^ = (w+ , ) . The V matrix was obtained as 



V h (9) 

V sm cos ^ ; 

where 6a is defined by 

cos 0a = , ^° ■ sin 0a = — , (10) 

Eigenvectors are independent of qa- They depend only on Eq. As varies from — oo through zero to cx), 9a varies 
from zero through 7r/2 to tt and varies from ^ g ^ through to ^]'^ • Thus, at = or at = Eq 

+ ^ cot A52 , two background eigenchannels are mixed equally. For this reason is regarded as the avoided 

crossing point energy. This energy does not coincide with the energy = qa where two eigenphase shift curves are 
separated most. Let U = U^V , then U diagonalizes the 5 matrix, that is, the transform USU is the diagonal matrix 
g~2i(5 rpj^g jj -fj^^^pjjj jg obtained from the V matrix by replacing 9a with 6*^ = 9a + 9'^ , where 0° parametrizes C/° 
matrix as 



t/°=r°^j , (11) 
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With the new parameters and Pauli's spin matrices, the S matrix was found to be expressible as 

S = e-'^^^'^+^-^'^O , (12) 

where 

< = zcos0; + £sin0;. (13) 

Smith's time delay matrix Q {=ihS^ ^ 

can be easily obtained by substituting Eq. ( p^ into its definition 

and was found to consist of three terms 

Q :^^{lTr+(T -h'^Ta + tT -h'jTf), (14) 

one due to the resonance, one due to the avoided crossing interaction, and one due to the change in frame transfor- 
mation as a function of energy, where 

n'^ = y X n'^ cos(5q - y sin(Sa, (15) 

and is orthogonal to n'^. 

The time delay due to the resonance takes a symmetric Lorentzian form 

^^dSJE) Ah 1 
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and the time delay due to the avoided crossing was found to take a form of a functional of the Beutler-Fano formula 

ra{E) = -ME)^ ^-'^ : (17) 

when..<g. ^^^^ 
where parameters r and are defined by 



Vr2 - Ar2 ^ ^ 

Ar ' (1^) 

g.= ^cotA°2, (20) 

The asymmetry of as a function of energy is brought about by the nonzero value of Qr which is proportional to 
the shift of the avoided crossing point energy from the resonance one. Thus the asymmetry of Tq is caused by the 
mismatch in the positions of the avoided crossing point and resonance energies. The time delay due to a change in 
frame transformation was found to take the following form fl^ ] 

rfiE) = rr{E)— (21) 

V/BF(er,gT) 

Because of the last term of Eq. (^4|) , eigenfunctions of the Q matrix are different from those of the S matrix. The 
eigentime delay sum, which is equal to Qu = TrQ, is obtained as 



Trg = Tr, (22) 



since Tr<T = 0. The consideration of the transforms tJ^SU^ and tJ^QU'^ , which will be denoted as S and Q, will be 
proved to be more convenient later for the geometric consideration. The transforms S and Q are the scattering and 
time delay matrices with the background eigenchannel wavefunctions as a basis instead of the asymptotic channel 
wavefunctions whose forms in terms of the new parameters and Pauli's spin matrices are the same as those of S 
and Q but with vectors fia and fif which are obtained from fi'^ and by replacing 9'^ with 9a- 

The connection of the time delay matri x Q with the time delay experienced by a wave packet was first considered 
by Eisenbud and extended by others ||l^. According to their work, Qu is the average time delay experienced by 
a wave packet injected in the ith channel. Here, the average time delays due to the avoided crossing interaction are 
given by cos 0^/2 and —Ta cos 0^/2. Similarly, the average time delays due to the change in frame transformation 
are — sin f?^ cos (5a/2 and sin0^ cos(5a/2. Time delays due to the avoided crossing interaction and the change in 
frame transformation are out of phase by 7r/2. Overall, Qu = ^(r^ + TqCos^^ — sin 0^ cos (5a) and Q22 = \{Tr - 

Ta COS 9a + Tf s'm 9'^ COS 5a ) ■ 

In analogy with the spin i system, the time delay matrix Q was expressed in terms of polarization vectors and the 
Pauli spin matrices as 

Q = ir,,(l + P, -cr + P^-cr), (23) 

where polarization vectors are defined by 

P„ = P, = ^n'/. (24) 

Like the spin i system, it was found that the absolute values of P^ and P/ are restricted to < |Pq| < 1 and 
< |P/| < 1. In the present case a complete depolarization means that eigentimes delays are the same regardless of 
eigenchannels, while a complete polarization means that eigentime delays are for one eigenchannel and Tr{E) for 
another eigenchannel. Eigenvectors for eigentime delays due to an avoided crossing interaction and due to a change 

in frame transformation are orthogonal to each other and contribute to the total eigentime delays as w + r| = 



7"rV|PaP + It "^as found that 

\Pa\' + \Pf\'^l. (25) 



4 



Since and P/ are mutually orthogonal and |Pap + |P/p = 1, we can define a vector Pt = Pa + P/, whose 
magnitude is unity. Its formula may be obtained straightforwardly but hardly used. Instead, the formula of its 
transform Vt — U^'PtU^ is exclusively used, which is much simpler and given as 




The transform Vt is the total polarization vector for the time delay matrix Q {—U'^QU^) with background eigen- 
channels used as the basis. The similar transforms Va and "P/ of Pa and P/ will be used later and satisfy the same 
relations Vt = Va + 'P f and 

|PaP + |P/f = l. (27) 

With the total polarization vector, the time delay matrix becomes 

Q= ir,(l + Pfa). (28) 

Since (Pt • cr)^ — 1, eigenvalues of Q or total eigentime delays are obtained as zero and r^, the time delayed by the 
resonance state. Though time delays due to an avoided crossing interaction and a change in frame transformation are 
asymmetric with respect to the resonance energy and therefore the energies of the longest lifetimes are not matched 
with the resonance energy, the energy of the longest overall eigentimes delayed is exactly matched with the resonance 
energy. 



III. PREPARATION FOR THE GEOMETRICAL REALIZATION 



In the previous work, some of the interesting things were noticed but could not be explained. Some of them are 
summarized below. 

• Why are eigenvectors of the S matrix independent of qa while its eigenphase shifts are not? 

• Why do the energy behaviors of 5a{E), Ta{E), and Tf{E) follow Beutler-Fano formulas? 

• Why does Ta take the Beutler-Fano formula in the energy scale of instead of Ca in contrast to the case of Sa 
though the former is obtained as the derivative of the latter. 

• Why is |PaP + |P/P = 1 satisfied? 

• What is the meaning of the parameter r^? 

In the previous work, we got some help by making an analogy of the system with a spin model, especially in interpreting 
the time delay matrix Q with polarization vectors Pa and P/ which are borrowed from the spin model. But the 
analogy with the spin model is not fully exploited. Here we show that by exploiting the analogy further, we can 
give the explanations of the above questions. In particular, we succeeded in giving the geometrical realization of the 
Beutler-Fano formulas. 

Before starting the geometrical realization of the previous results, let us first rewrite some of the previous results 
suitable for the geometrical realization. 

First, we notice that Eqs. (||) and ( p7| ) are simpler than the corresponding Eqs. (0) and (|l8|). This indicates that 
the square root of the Beutler-Fano formula (^ seems to be more fundamental than the original one. Next we notice 
that Eq. (^) resembles cot 6r = — and cot^a = —^a- Thus the square root of the Beutler-Fano formula may be 
regarded as an energy parameter eBF(e, 9, ^°)- Then Eq. (^) takes the suggestive form 

cotJa = -eBF(ea,9a,A°2), (29) 

where 

eBF(6a,<?a, A;^) = COt A^^ ^f=J% (30) 
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{AI2 is the value of 6a at Ca —00). But there is a one drawback when the square root of the Beutler-Fano formula 
is considered as an energy parameter. It is not a monotonically increasing function of energy. It has a minimum when 
q > and a maximum when q < 0. Hence ebf will be considered here merely as a convenient notation. 

Eq. (^ ) suggests another angle 9f satisfying Va — na cos9f and V f — hf sm9f. Its cotangent is obtained as 



cot Of — — 



1 er -qr 



(31) 



Eq. ( |3l| ) indicates that l/r becomes cot Of at Cr — > —00. The angle of Of at — *■ —00 is identified with the angle 
which the polarization vector Vt or fit makes with fia- That angle will be denoted as Ot- Eq. ( |2^ ) shows that the 
angle Ot is obtained as 



cos t)t 
sin Ot 



AT 

Vr2 - Ar2 



(32) 



and with it the spherical polar coordinate of fit is given by {I,0t,—Ai2 



cot Of 



' cott 



Now with 6t, Eq. (31) becomes 
eBF(er, qr,Ot). 



(33) 



With the new angle 9f, Ta becomes r^cos^/, which explains the complicated form of as a functional of the 
Beutler-Fano function in contrast to that of cot 6a- 
As a result of rewriting, we obtain four equations 

cot 6r — — Er, 

cot 6*0 = -ea, 

cot 6a = -eBF(ea,ga, A?2), 

cot 6*/ = -eBF(er, gT,6't). 



(34) 



The use of the geometrical parameters, 6r and Oa, in place of er and makes the geometrical realization of dynamic 
relations possible. Our aim is to obtain the dynamic formulas for SBpi^a, qa, A12) and eBF(er, 9r, ^t) by converting 
the geometric relations containing 6a and 9f back into dynamic ones. We will sometimes abbreviate eBF(fa5 ^aj ^12) 
as eBF,a and eBF{^r,qT,Ot) as eep.r- Before ending this section, let us note the following formulas for the line profile 
indices qa and q^ 



qr = 



cot 6a{ea = 0) 
cot 6a{ea -00)' 

cot 61/ (e^ = 0) 



cot Of{er 



-00) 



(35) 



They can also be expressed as 



cot 5a , TT 

cot A\2 2 

COt^f TT 

qr = — — - when Or = -• 
cot Ot 2 



(36) 



IV. GEOMETRICAL REALIZATION 



The geometrical realization is based on that to each unimodular unitary matrix in the complex two-dimensional 
space, there is an associated real orthogonal matrix representing a rotation in real three-dimensional space ||l^. The 
general two-dimensional unimodular unitary matrix can be written as e~'2'^ " as its determinant can be easily shown 
to be unity using Tr(cr) = and thus unimodular by the definition of unimodularity. The associated real orthogonal 
matrix will be denoted as Rn{0), the rotation matrix about the vector fi by an angle defined in an active sense. 
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Let us first consider tiie S matrix. It is unitary but not unimodular [det(S') 7^ 1] and can not be associated with 
the pure rotation alone. But after extracting det(S') which is equal to e""'^ for isolated resonances (a similar formula 
holds for overlapping resonances, where Sr is replaced by the sum over ones from all resonances [|9|), the remaining 
part of the scattering matrix will be unimodular and may be associated with a pure rotation. According to Eq. (|l^), 
the remaining part is e^"'"'^ "" and may be associated with the rotation about the vector fi^ by an angle 2Sa- We will 
explore the possibility of this explanation of the S matrix in below by deriving Eq. ( [l^ ) in a more systematic way. 

In the previous section, the S matrix was diagonalized by two unitary matrices and V given by Eqs. ( pi] ) and 
(^. Actually the theorem in Ref. limits that t/" and V are real orthogonal since the S matrix is symmetric. With 
them, it is rewritten as 

S = U°Ve~'^'^VU°. (37) 
Note that the two unitary transformations U'^ and V can be written in terms of Pauli spin matrices as 

V^e-'^^-y. (38) 

Notice that argument matrices of two exponential functions are commute and therefore U'^V — e-^ki^a+s°)"'-y ^ p^^ 
before, let us denote 9a + 6*° as 9'^. The diagonalized matrix e~^*^ of the S matrix can be expressed in terms of Pauli 
matrices as 



e-'''=[ " )=e-(^-^+^-^). (39) 



Substituting Eqs. ( |38| ) and (39) into Eq. (|37D, we obtain 

g — g-j%-CT-l(g-i(<5sl+<5aCr-2)gi%-<T-y^ J-^Q-j 



In order to give the geometrical interpretation to Eq. (40), a long preliminary exposition is necessary. Let us start 



with considering cr ■ r and transform it into a new matrix cr • r' by a general 2x2 unitary transformation e *2'^ "- as 
follows 

e-4'^-"cr.re'^'^-" = (tt'. (41) 
The left hand side of Eq. (^) can be calculated using ||2l[] 

g.Sog-.5 ^ o + 6] + ^[^, [5, 6]] + ^[5, [5, [S, 6]]] + . . . (42) 

and (cr • a)(cr • b) = a • b + icr • (a x b) and the result is that r' is just the vector obtained from r by the three 
dimensional rotation matrix Rn{9) about the vector n by in an active sense as 

v'^Rf,{9)Y. (43) 

Only in the form of the similarity transformation (^) the homomorphism that a 2 x 2 unimodular unitary matrix 
is associated with a three dimensional rotation holds. According to this interpretation, the unitary transformations 
U'^ and V for the symmetric S matrix ( |37| ) correspond to the three-dimensional rotations about the y axis through 
angles 9''^ and 9a, respectively, and their overall effect U'^V is equal to the rotation about the y axis by 9'^^ = 9a + 9''\ 
Therefore, the original frame transformation for the symmetric S matrix becomes as a rotation about the y axis in 
the "hypothetical" real three-dimensional space. This hypothetical real three-dimensional space is different from the 
Hilbert space and called the Liouville space. It is the space spanned by the set of vectors Gx, Ty, and az which are 
orthogonal in the sense that 

Tr(a,aj) =25„-, (44) 

and extensively studied in Ref. [ p2[ . Any traceless 2x2 Hermitian matrices, for example h, can be expanded in 
this vector space as h — xa^ + yay + za^ = {x, y, z) = a ■ r. We can lift the restriction of traceless in Hermitian 
matrices if we include the unit matrix 1 as another basic vector in addition to ax, cry, az- Then the three-dimensional 
Liouville space is a subspace of this four-dimensional Liouville space. Note that two subspace {1} and {ax,ay,az} are 
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orthogonal and therefore either the subspace {a-x, cry, cr^} or {1}, or the whole space may be chosen freely depending 
on the situation without having any trouble. 

Now, Eq. (|4^) can be viewed in two ways. It can be viewed as a rotation of the vector r into r' by the rotation 
matrix Rn{9) as expressed in Eq. (^). Or it can be viewed as the transformation from the xyz coordinate system 
to the x'y'z' coordinate system by the rotation matrix Rn{—0). The latter view, though obvious, can be shown to 
be true using the following mathematical transformation. Let us regard cr and r as the column vectors. Then the 
scalar product cr • r' can be written as a matrix multiplication of the row vector cr with the column vector r', namely, 
cr • r' = or'. The support for the view of the coordinate transformation is obtained by the following transformation 

o-r' = aRn{e)r = [R^)a]r. (45) 

Since the diagonalization of the operator cr • r yields its eigenchannels, the vector r in the three-dimensional Liouville 
space is enough to uniquely specify the eigenchannels of the traceless Hermitian matrix cr-r. Conversely, eigenchannels 
may be regarded as a vector in the three-dimensional Liouville space. 

Since any real orthogonal frame transformation is of the form like Eq. (^), it may be generally said that a real 
orthogonal frame transformation in the complex two-dimensional space corresponds to a rotation about the y axis 
in the real three-dimensional Liouville space. Since the matrix corresponding to any 2x2 Hermitian operator is 
diagonal in the basis of eigenchannels by definition of eigenchannels and can be written as al + bcr ■ z, a dynamical 
process along an eigenchannel corresponds to a process along the z axis in the real three-dimensional Liouville space 
and leads to a variation in length of the vector. Thus the y axis in the real three-dimensional Liouville space can 
be regarded as the axis for the real orthogonal frame transformations and the z axis as the axis for the dynamical 
processes along eigenchannels. 

We have theorem that Ce^C~^ = e^^^ which can be easily proved by using Eq. (|4^) . Using this theorem and 
Eq. (^), we have 

Using Eq. (^) and Eq. (Hj), Eq. becomes Eq. with fi'^ now interpreted as 

h',^Ry{e',)z. (47) 

Or nj, can be regarded as the z' axis in the x'y'z' coordinate system, i.e., fi'^ = z' . Let S = e^'^*^ , i.e., A' — 
\{5t. 1 + 5a(T ■ njj). A' is a vector in the four-dimensional Liouville space. Or, if we exclude the isotropic part in A', 
^(5acr • nj, is a vector in the three-dimensional Liouville space which may be obtained by rotating the vector ^Sacr ■ z 
about the y axis by an ang le 9'^. The \8y,1 term in A gives the phase shift owing to the isotropic influence of the 
background potential scattering and the resonance on eigenchannels. Likewise, the ^(5aCr • term gives the phase 
shifts owing to the anisotropic influence of the background scattering potential and the resonance on eigenchannels. 
Therefore, the length ^5 a of the anisotropic term in the three-dimensional Liouville space denotes the degree of 
anisotropic influence on eigenphase shifts by the background scattering and the resonance. 

Let us now consider the time delay matrix. If the time delay matrix is written in Lippmann's suggestive form 
|p3| , Q = S'^tS (t is the time operator defined by ih-^), it is apparent that the unitary matrix which gives the 
similarity transformation is now the matrix and can be associated with the rotation matrix according to the 
theorem mentioned above when det(S') is extracted from it. Using the relation 



dz 



- (e^^) = A^'e^' = e'^M^ (48) 



Eq. (^ ) is easily differentiated with respect to energy to yield 



By multiplying with the adjoint of ([40[), the Q matrix becomes 

where use is made of the fact that S is unitary and thus S^S ~ 1 for the first and third terms. The matrix 
multiplication in the second term of the right hand side of Eq. ( pO| ) is just o ■ n'^ as was already done before. The 
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sum of the first and second term is the time delay due to the energy derivatives of the eigenphase shifts and called 
"the partial delay times" by some group of people j|] . 

The parenthesized part of the third term is the time delay due to the change in frame transformation and has an 
interference effect between two contributions, one due to the change in frame transformation from the asymptotic 
channels to the background eigenchannels (iPe^'^^I^'e) (^'^ ~ lj2,...,n) and the other due to the change in frame 
transformation from the background eigenchannels to the asymptotic eigenchannels (V'b'' I V'_b *■'"'') (^w — 1,2,. ..,n). 
The change in frame transformation does not take place in the direction of the rotation axis given by y but in the 
rotation angle since the rotation y axis is fixed in the Liouville space as energy varies. The first contribution has the 
term S~^(t ■ yS which is the similarity transformation of the operator a ■ y by S. Substituting Eq. (^ ) for S, this 
term becomes 

S+CT ■ yS = e-'^'^-Se*'^"'^-^ • ye'^'^-'^-^e*^'^-*, (51) 

where use is made of that 2 y and cr • y are commutative. The scalar factor e in S does not appear in Eq. 
( ^ ) as it is multiplied by its complex conjugation in to become unity. According to the theorem, the unitary 
transformations in the right hand side of Eq. ( ^l| ) correspond to two consecutive rotations, at first about the z axis 
by —2da and then about the y axis by Oa- By the first rotation, cr • y becomes cr • (xsin25Q + ycos2Sa)-, i.e., 

e^s.cT-z^ . ye-'*"-^-^ =T-[R, i-2Sa)y] = cr ■ {x sin 2Sa + y cos 2Sa) ■ (52) 

By substituting Eq. (63) into Eq. (pW) and cr • y's being replaced with e~*~'^'^cr • ye^~'^'y , 

S+cT ■yS-(T-y = e-' "-^''-y [a ■ [x sin 25a + y cos 25a) - (t ■ y] e*^'"'^ 

= e-'^'^-y [2 sin 5a(T ■ {x cos 5a - ysin5a)] e'^""-*. (53) 

The bracketed part of Eq. ( |5^ ) is equal to the rotation of the x axis about the z axis by —5a multiplied by 2 sin 5a, 
which is the overall effect of the interference. Fig. |l] shows this process of interference as a vector addition in the 
three-dimensional Liouville space. The time delay due to the change in frame transformation, the third term of the 



right hand side of Eq. ( |50[ ) , becomes 



hsin5a^e~'-'' y [x cos 5a -ysin5a) e'-'^'^ 

= hsin5a^e-''-^'' ye'^''-^(T ■ ie^^^'^-^e*^'^-^ 
dE 

de' .6a ; , Sa 

= hsin5a-r-e'-'' ' a- ■ x' 

dE 

= hsin5a'^(T-x", (54) 
dE 



where the second equality is obtained by applying Eq. ( |4l|) tw ice to the matrix term to obtain cr • [Ry{9'a)Rz{—5a)x 



which becomes cr • Rz'{—6a)x' and then by applying Eq. again. In the last equality of Eq. (54), we introduced 
another new x"y" z" coordinate system which is obtained from the x'y'z' coordinate system by tl: 
in the passive sense. In the active sense, x" = Rz'{—5a)x' . 

Substituting Eq. (m) into Eq. (pOh, the time delay matrix Q is obtained as 



le rotation Rz'{5a) 



Q = n[l—+a.z — + cr.x sm<5.— j, (55) 

and is equal to Eq. ( [l^ ) when z' and x" are identified with the unit vectors h'^ and h'^ (n^/^ and fi^, in Ref. 0) and 

hdSs/dE {— hd5r/dE), hcr-z'd5a/dE, and h sin 5a(T-x"d9a/dE are identified with ^r^, ^t^Pq and ir^P/, respectively. 
By substituting Pa = z' cos Of — z" cos Of and P/ = a;"sin6'/, Eq. ( |55| ) can also be transformed as follows 

Q = ^Tr[l + (T ■ {z" cosOf + x"sin9f)] 

1 / .// -Of ~, 

= —Tr ll + e 2 "cr-ze^ " 
= ir,(l + ^.z"'), (56) 
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where still another x"'y"' z'" coordinate system is introduced. In the active sense, z'" = Ryii{Of)z" . Eqs. ( JSq ) and ( |5q ) 
tells us that time delay matrices due to the avoided crossing interaction, due to the change in frame transformations and 
the total time delay matrix take simplest form in the x'y'z', x"y"z" , and x"'y"'z"' coordinate systems, respectively. 

Eq. js^ ) equals Eq. (|28| ) and therefore z'" equals Pt- The vector z"\ and accordingly Pt, can be obtained from z 
by successive rotations by 

P, = z'" = Ry„{0j)R,,{^da)Ry{e'Jz = Ry„{ef)Ry{e'Jz. (57) 

As mentioned before, it is better to consider Vt = U^'PtU'^ rather than P( itself since the formula for the former is 
simpler that that for the latter. Vt is the polarization vector pertaining to Q = U^QU^ which is the time delay matrix 
in the basis of background eigenchannel wavefunctions. This suggests that it may be better to take the background 
eigenchannel wavefunctions rather than the asymptotic channel wavefunctions as a starting channel basis. From now 
on, let us redefine the xyz coordinate system as the coordinate system pertaining to the background eigenchanncls. 
Let us use the name of x^y^ z'"'' coordinate system as that pertaining to the asymptotic channels. Definitions of other 
coordinate systems remain unchanged. With this redefinition of notation, the formulas for S and Q corresponding to 
Eqs. ( |l^ ) and (14) becomes 



g _ g-j%-CT-t/g-i(5El+5„CT-i)gi%-o-v „ g~i(i5sl+(5„<T-n„)^ j-gg-j 
Q^]^{Trl + Ta(T -fla+TfCT -hf) = {1 +V a ■ (T +V f (t) = ^T^ (1 + • Ct) , (59) 



with 



ha — z COS 9a + X sin 6a = Ry{9a)z, 

hf :^ y X haCOsSa - {jSmSa = (-(5a)x'. (60) 

In place of Eq. (|5^), we have 

Vt^ht^ Ryn{0f)R,,{-6a)Ry{ea)z = Ryn{Of)Ry{ea)z. (61) 



By substituting the relations cot^a = — ea and cot Of = — CBF.r into Eq. (|61|), it is checked that the same formula as 
Eq. (|2^ ) is obtained for fit- Note that the formula ( P6| ) for fit is independent of energy in contrast to ha (hf) which 
varies from z (x) through x {—z) to —z {—x) as energy varies from — oo to oo. This holds generally at least for the 
multichannel system in the neighborhood of an isolated resonance and derives from the fact that only one type of 
continua can interact with a discrete state (see Eq. ( |A6D in Appendix ^ and Rcf. jl^ for more general systems). 

So far, several coordinate systems are considered such as xyz, x'y'z', x"y"z", and x"'y"'z"' coordinate systems 
pertaining to the eigenchanncls of 5°, 5' or cr • Va, ^ ' ^ a-nd a ■ Vt, respectively. These coordinate systems are 
shown graphically in Fig. ||. According to Eq. (|3^), the spherical polar coordinate of ht is given by {l,Ot, —^12) in 
the xyz coordinate system. Since ha lies on the zx plane, the absolute magnitude of the azimuth of ht is equal 
to the dihedral angle between two planes whose normals are given by z x and z x ht- Let us now consider the 
coordinate of ht in the x'y'z' coordinates, where z' = ha- The angle which ht makes with z' is 0/ and the azimuth of 
ht may be obtained by considering the z'x" {=z"x") plane. Note that x" is equal to hf and ht lies on the z'x" plane 
meaning that the azimuth of ht is identical with the dihedral angle which the z'x" plane makes with the z'x' plane. 
Since x"y" z" coordinate system is obtained from x'y'z' coordinate system by rotating about the z' axis by —6a, the 
dihedral angle which the plane z'x" makes with the z'x' plane is 6a- Therefore, the spherical polar coordinate of ht in 
the x'y'z' coordinate system is (1, 9f, —6a)- See Fig. ^ to understand the explanation graphically. Since the dihedral 
angle between two planes whose normals are ha x x' and ha x ht is 6a, the dihedral angle between two planes whose 
normals are given by x z and ria x is tt — (5^- With this, we can construct a spherical triangle AAPQ with 
vertices formed with the endpoints of z, ha, and ht, where the vertex angles opposite to the edge angles 9f and 9t are 
A^2 ^'Hd TT ~ 6a, respectively, as shown in Fig. ^ The vertex angle opposite to the edge angle 9a can be shown to be 
Sr by making use of the following relation 

^-iAl^rr-z ^-iSr<T-nt _ ^-iSarr-fba ^g2) 

which is the spin model version of the relation S = S'^{tti, + e^^^'^'^TTa) (see Appendix |^ for the derivation). According 
to Appendix Eq. ( |6^ ) can be expressed using the rotation matrices in the Liouville space as 

R,{2A°2)RAA2Sr) = RnA'^Sa) (63) 



10 



and manifests that the vertex angle opposite to 6a is Sr- 

The dual spherical triangle of AAPQ may be constructed by converting vertices of the original triangle to its edges 
and edges of the original triangle to its vertices, according to the rule described in Ref. Q. According to the rule, the 
vertex angles of the dual spherical triangle are obtained from the corresponding edge angles of the original triangle by 
subtracting the edge angles from vr like n — 9a, it — Of and n — 0t and the corresponding edge angles are also obtained 
similarly like tt — Sr, tt — and Sa- The dual spherical triangle constructed in this way is shown in Fig. |^. 

Before considering dynamic aspects of the laws holding for the spherical triangle, let us comment on Eq. (|6^, or 
the equivalent 5°(7rb + e'^'^^'''Ka). For this purpose, let us define phase shift matrices A*', A^, A by 5° = e~^^^ , 
TTh + e~^^^''na = e^^*'^'^, S = e^^''^. Phase shift matrices are easily obtained as 

A0 = i(5^1 + A;2<T.z), (64) 

Ar = ^{Srl+SreT -fit), (65) 

A = ^{5^1 + 6a<T-na). (66) 

For Eq. (^), if — 0, two eigenchannels have the identical background eigenphase shifts = 62 — The 
eigenchannels are isotropic with respect to the potential that brings about the background phase shifts. When two 
eigenchannels react differently or anisotropically with respect to the potential, eigenchannels have different phase shifts 
61 S2, or AI2 0. (The off-diagonal term of 5*° gives the transition amplitude and is thus caused by the channel- 
channel coupling. Since the off-diagonal term is implicitly included in eigenchannels, the potential that eigenchannels 
feels includes the channel-channel coupling effect.) The anisotropic term of (Q) contains the information on this 
phase difference and the eigenchannels. 

The phase shift matrix ( |6^ ) is a vector whose coordinate is ((5^, 0, 0, Aj'j) in the four-dimensional Liouville space. 
Or, if we consider only the anisotropic term, it is a vector in the three-dimensional Liouville space, whose magnitude is 
■iA52 and whose direction is z. Though background and resonance scattering contributions appear as a single product 
term in Eq. ( |6^ ) in the S matrix, two contributions are not simply combined in the case of the phase shift matrix A. 
For the isotropic parts of two contributions to A, the combining rule is simple and they are simply added up to give 
the isotropic part of the phase shift matrix A as — ^{S^ + Sr)- The combining rule of anisotropic terms is not 
so simple. According to the Campbell-Baker-Hausdorff formula the anisotropic part of the phase shift matrix A 
is expressed as a very complicated infinite sum of multiple commutators of the anisotropic parts of A'^ and A^ as 

2A — S'sl = A^jC • z + 5r<y ■ fit — — [A"2C • z, Srcr ■ fit] 

- ^ ([[A"2cr • Z,5r(T ■ ht],SrCr ■ nt] - [[A"2cr • Z,5r(T ■ nt], A"2cr ■ z]) -\ (67) 

But, the geometrical construction in the Liouville space provides a simple combining rule. For that, at first we ignore 
the magnitudes of vectors in the three-dimensional Liouville space and consider the spherical triangle made up of 
endpoints of the unit vectors corresponding to the anisotropic terms. The magnitudes of vectors corresponding to 
anisotropic terms, instead, are utilized as the edge angles of the spherical triangle. This is the procedure we take 
when we interpret Eq. ( |62| ) as giving the remaining edge angle Sr- Trigonometric laws of the spherical triangle then 
provide a details of the combining rule which are the subject of the next section. 

Let us comment one more thing on the avoided crossing interaction. If we use background eigenchannels as a 
basis, then the 5° matrix is already diagonal by definition in that basis. Then for processes occurring along this 
eigenchannels, there is no channel-channel coupling. If a discrete state is included into the system, background 
eigenchannels are no longer decoupled and interact each other through the indirect continuum-continuum, or channel- 
channel, coupling via the discrete state. This indirect channel-channel coupling brings about the avoided crossing 
interaction in the curves of eigenphase shifts of the S matrix. Therefore, the avoided crossing interaction is not 
devoid of the resonance contribution. What is devoid of in the avoided crossing interaction is the isotropic resonant 
contribution. It includes the anisotropic resonant contribution. 

V. CONNECTION OF THE GEOMETRICAL RELATION WITH DYNAMICS 

Now let us describe the dynamical aspects of the geometrical laws holding for the spherical triangle, such as the 
laws of sines, the laws of cosines. Cotangent laws including four successive parts and laws including five successive 
parts (see Ref. p^) are derivable from the laws of sines and cosines but deserve a treatment as separate laws. 
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Let us consider the cotangent law 

sin cot 6a — — sin 6a cot 6t + cos 9a cos A°2 . (68) 

When cos 6*0 — —ea/ + 1 and sin 6*0 = l/-\/e^ + 1 are inserted (for the sign convection, see [^), Eq. (68) can be 
put into Beutler-Fano's formula (|^) for cot Sa as follows 

cot Sa = . ''■0 
sm 




= -cotA;2^^=J^, (69) 

where Qa is identified with — cot 6t / cos A!|'2 and can be easily checked to be equal to the previous definition (^ with 
the use of Eq. (|2|). 

One of the cotangent law for its dual spherical triangle is given by 

cot 9f sin 9t — — cos(7r — Sr) cos 6t + sin(7r — Sr) cot Ai2- (70) 

With cos 5,. = — e,./y^e2"+T and sinJ^ = l/V^r + 1j Eq. ( [70| ) can be put into Beutler-Fano's formula ( ^ for cot6'/ 
as follows 



cot 9f — -r— TT (cos Sr COS 6*4 + sin (5^ cot A?n) 
smt't ^ ' 

I ( er 1 , n ^ 

— ;r — cos 9t cot AVo 



where qr is identified with cot Aj'2/ cos^f, again equal to Eq. (20). 
Using Eq. ( |3^ ) and the convention ||2^, the sine laws, 



= -cot0,i^, (71) 
\/H + 1 



are translated into 



sin A'j'2 sin Sr 

sa\9f sin 

sin Sr sin Sa 

sin 9a sin 9t 



(72) 



e2 + l = sin2A?2(e^ + l)(e|F,r + l), (73) 
e2 + l = sin2et(e2 + l)(4 +1), (74) 



respectively. Eq. (|7^) was used to derive Eq. ( |2l| ) in Ref. [Q. 

(Equating the a;, y, z components of ht with those of fia cos 0/ + n/ sin 0/ yields the law involving 5 successive parts 
sin6'tCosA52 = sin 6*0, cos 6*/ + cos 6*0 sin 6*/ cos (5a , one of the laws of sines, sin A5'2/sin6'/ = sin (5a/ sin ^t, and one of 
the laws of cosines, cos9t = cos 0a cos0/ — sin 0a sin0y cos Sa, in that order. Though such an equality looks irrelevant 
to the laws of the spherical triangle at a first glance, it actually has to do with the laws of the spherical triangle since 
it is used to obtain the vertex angle tt — (5a of the endpoint of fia-) 

The presence of the dual spherical triangle indicates that there is a symmetry with respect to the exchange of 
vertices and edges. The comparison of the spherical triangle in Fig. ^ with its dual one in Fig. || shows that the 
exchange of Sr, 9f, and Aj'2 with tt — 0a, <5a, and tt — 0t transforms the spherical triangle into its dual and vice versa and 
thus any trigonometric laws will be invariant under this exchange. Besides the geometrical laws, other laws containing 
not only geometric parameters but also other types of parameters should remain as valid expressions with respect to 
this exchange. With this requirement, in order for cot Sr = —tr to remain as a valid expression with respect to the 
exchange, should be replaced by — ea- The right hand side of qr = cot A]'2/ cos0t becomes cot 9t/ COSA52 which 
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is equal to — Thus <?r is replaced by —qa under the exchange. Similar procedure shows that esF.a is transformed 
into eBF,r under the exchange. The variables with their conjugated ones are summarized in Table [[ This symmetry 
under the exchange yields many relations without derivation thus saves a lot of efforts. It can also be used to check 
the validity of derived equations. Let us take a few examples. If it holds the relation that 



j_ 

Qa 



sin dt sin A°2 

another valid relation is obtained as 



sin sin 9t 

by exchanging variables according to Table ^ If it holds that 



(75) 



(76) 



then the following 

-cos (5a (78) 



can be obtained by the same procedure. 

Geometrical realization reveals that complicated behaviors of dynamical parameters like 5a and Of as a function 
of energy are nothing but the result of a simple geometrical traversal along the great circle shown in Fig. Before 
examining the behaviors of dynamic parameters as functions of energy, it is noted that the vectors z and fit are fixed 
in the real three-dimensional Liouvillc space while fia changes its direction as energy varies. The constancy of the z 
vector derives from the usual assumption of energy insensitiveness of background scattering. The cons tanc y of the 



fit vector derives from that the time delay matrix has only a resonant contribution as shown in Eq. (A6) and its 

eigenchannels consist of Fano's energy independent ij)^^^ and continua orthogonal to it. As the energy tr varies from 
— oo to cx), 6a undergoes a change from to tt which corresponds to the semicircular traversal of the point P from A 
to the opposite point —A along the greatest circle while the points A and Q keeps fixed in Fig. |[ In this semicircular 
traversal, the angle tt — 5a varies from tt — A°2 when P coincides with A to when P coincides with —A. The angle 
5a, accordingly, varies from A?2 to TT - A?2. (The angle Of varies similarly from 9t to n — 9t.) Traversal enjoys a 
special symmetry when 9t = 7r/2. Let the point A be taken as a polar point. Then the side PQ becomes part of 
the equator at the middle of the traversal where 9a becomes 7r/2. Since any meridian makes the right angle with the 
equator, the angle 5a which the chord AP makes with the equator becomes a right angle, i.e., 5a = Tr/2. Now let us 
consider the deviation of the point P from the equator. Let 5a = 7r/2 + y at 0a = 7''/2 + x. Then, by the symmetry 
of the spherical triangle, 5a — tt /2 — y at 9a ~ t: /2 — x. Napier's rule 

cot 5a = cos 9a cot A'j'a (79) 



holding for the right spherical triangle satisfies such a symmetry |g8|. Obviously, qa — ^ for Eq. (|79|). When 9t ^ 7r/2, 
5a is no longer 7r/2 when 9a — 7r/2. The occurrence of the mismatch in energies where values of 9a and 5a become 
7r/2 amounts to the addition of sin^Q term on the right hand side of Eq. ([tq]), which causes the value of qa to deviate 
from zero. The value of qa can be obtained as — cot 9t/ cos A^j by substituting — cot 9t/ sin A52, one of Napier's rules 
holding when 9a = 7r/2, for cot 5 a into Eq. (^6|). It can be roughly stated that the asymmetry of the Beutler-Fano 
formula for cot 5a derives from the asymmetry of the geometry. 

The concurrent change oi 9f with the increase of the arc length 9a as P traverses can be obtained by differentiating 
the cosine law cos9f — cos 9a cos 9t + sin 9a sin 6*4 cos Aj'2 with respect to 9a keeping 9t and fixed, which becomes 

d9 

— sm9f — p — — sin 9a cos 9t + cos 9a sin 9t cos ■ (80) 
dOa 

The right hand side of Eq. ( |80| ) becomes — sin0/cos(7r — 5a) according to the law containing five successive parts, 
which finally yields 

^=-COS^a. (81) 
d9a 
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Similar derivatives are obtained as 



dSa „ 

dSa ,/) • r 

— - = cot&f smda, 

dcot 0a ^ sin 
d cot Sr sin 9t 

In the spherical triangle AABC, Gauss-Bonnet theorem becomes 

area(AABC) 

which states that the sum of interior angles of a spherical triangle exceeds tt by the solid angle Q defined by 
area{AABC) / . In the present case, the sum of interior (vertex) angles is tt + 6r + — 5a- Hence, the solid angle 
O is 

f7 = + A°2 - (5a = 5r + (5s - 2(5§ - (5a = 2(5_ - 5l). (84) 

The solid angle of the spherical triangle AAPQ is easily calculated as 2A52 when the point P coincides with the 
antipode of A. Then the solid angle of the spherical triangle varies from zero to 2A]'2 as the point P varies from the 
point A to the point —A and, accordingly, 5- varies from ^2 to (5j as energy varies from —00 to 00, which is consistent 
with the result of Ref. . 

So far, a geometric realization of the S matrix and Q matrix has been considered. Let us now go back to the 



original questions we had in the beginning of Sec. Ill and see whether we can explain them. 

The first of the questions was that the energy variations of the eigenvectors of the S matrix are independent of 
qa while those of its eigenphase shifts, or more specifically (5q, depend on it. Let us start with that the eigenvectors 
or eigenchannels of the S matrix are obtained by the frame transformation of the background eigenchannels. Since 
the background eigenchannels are fixed, the energy variations of eigenchannels of the S matrix completely come from 
the frame transformation which in this case is parametrized with Oa- In Fig. H, 9a is the edge angle opposite to the 
vertex angle 5r- If Cr varies, 5r varies according to cot (5^ = — Cr. Since the edge angle 6a is the opposite to i5r, 6a may 
be expected to vary linearly with 5r- Such an expectation does not come out right. Instead, cot 0a varies linearly 
with cot 5r according to one of the relations in (|8^), dcot 6a/ dcot 5r = sin A'^2/sin0t, which is fixed in energy. The 
relation tells us that cot^a has a linear relation with but cot 0a may not be zero when = 0, in general. But, 
we can always introduce a new energy scale, let us call it ta, where cot 0a is zero at = and the proportionality 
constant can be set so that cot 0a = —Ca- The argument proves that cot 0a does not need no further parameter like 
qa- The reason why the energy variation of 5a needs the line profile index is already considered around Eq. ( [79| ) and 
need not be repeated here. 

If the second question which asks why the energy behaviors of 5a, Ta and Tf follow the Beutler-Fano formulas is 
changed like "is it possible to show geometrically that their behaviors follow the Beutler-Fano formula?" , the answer 
is yes and their behavior is the result of the cotangent laws holding for the spherical triangle, as we have done in this 
section. 

Let us answer the third question why Tq takes the Beutler-Fano formula in the energy scale of instead of ea 
though Ta is obtained as the derivative of ea- Note that the question on Tq can be paraphrased to that on cot 0/ since 
Ta — Tr{er){l + cot0^)"^/^. There are two such cotangent laws for cot0/ as follows 

cot0f — — ( cos 5r cos 6t -I- sin(5,. cot A?2) , (85) 
sm0t ^ ' 

cot 0/ = ~^~7r (~ cos (5a cos 0a + SluSa COt A'j'2) . (86) 
Sm (7a 

Eq. ( ^sl) expresses cot 0/ in terms of €r while Eq. (|86| ) expresses it in terms of €a- In Eq. (^5|), 5r is the only parameters 
which is a function of energy while, in Eq. (|8^), not only 5a but also 0a are parameters which are functions of energy. 
Eq. (|5|) gives the Beutler-Fano formula as a function of as we already saw in Eq. ([7l|). Eq. ( ^6| ) might also give 
a Beutler-Fano formula as a function of ea if 0a were a constant of energy but it fails to do so as 0a is a function of 
energy, too. On the other hand, if 6t were a function of energy, Eq. (|8^) could not give the Beutler-Fano formula, too. 
This argument reveals that cotangent laws of a spherical triangle alone are not sufficient to guarantee the presence of 
Beutler-Fano formulas. 
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Let us consider the answer to the fourth questfon which asks the reason why |Pt| = |Pa| + = 1 satisfied. 

Eq. ( [A6| ) shows that non-zero resonant behavior of time delay occurs only when the system is in the ■0^'' state, which 

derives from that the i/;^'* state is the only type of continuum which can interact with the discrete state. On the 
other hand, the unit magnitude of the polarization Pt means that only one continuum shows a resonant behavior 
while others do not. Thus Eq. (A6) proves that |Pt| = 1. 

The answer to the fifth question is provided by the identification of with tan^j and need not be considered 
further. 



VI. APPLICATION TO THE TRIATOMIC VAN DER WAALS PREDISSOCIATION 

Ref. and the present paper have developed the theory for the behaviors of eigenphase shifts and time delays. Let 
us now consider the application of the theory to the vibrational predissociation of triatomic van der Waals molecules. 
The theory can be applied in two ways. When the S matrix is known as a function of energy either experimentally or 
by a theoretical calculation, eigenphase shifts can be calculated directly by its diagonalization. Similarly, eigentime 
delays can be calculated from the S matrix. For these data, formulas for eigenphase shifts and time delays derived 
from the theory can be used as models with parameters in the formulas viewed as adjustable ones. Extracting best 
values of the parameters may be tried by fitting the data of eigenphase shifts and eigentime delays obtained by the 
diagonalization of the S matrix to the theoretical models. On the other hand, by using the formulas for the parameters 
themselves derived from the theory, parameters can be directly calculated without doing data fitting. Parameters 
obtained in two different ways, namely, by data- fitting and by using the theoretical formulas are not identical as 
the theory developed so far relies on the assumption that background eigenphase shifts and partial decay widths are 
constants of energy, which is usually a good approximation but does not hold exactly in the actual system. 

The data-fitting will be done only for the eigentime delay sum (^2|) and partial delay times 2hd5+ / dE and 2hd6-/dE. 
Eigenphase shifts will not be used for the data- fitting since they need the information on Ea and Fq, which is not 
available before the data-fitting. The data-fitting of partial delay times to the theoretical formulas 

^t,d6± .fdSr,d6a\ 1. , X 1 n , a \ 

2h = n ± ~ -(Tr ±Ta) = -Tril ± COSt'f ) 

dE \dE dE 2^ ' 2 ^ ^' 



2h 1 



It 



V(er--Zr)2+tan2 Otil + el) 



(87) 



can be easily done on the other hand since the information on Eq and F which are necessary to convert E to needed 
for Eq. ( ^7|) can be easily obtained from the data-fitting of the eigentime delay sum (|2^). (Eigenphase sum can also 
be used to get Eq and F.) 

Graphs of partial delay times are shown in Fig. |^ for several values of line profile indices. Some general characteristics 
of the graphs are noticed. 

1. Graphs of partial delay times 2hdS+/dE and 2hdd-/dE meet at e,. = Qt- 

2. As \qr \ cot A5'2, 



cos Of{er) 



1 when er < Qt 
-1 when er > qr, 



(88) 



and the partial delay times becomes 



2h 



2h 



dd+ 

Ile 

d5_ 

' dE 



Tr{er) when < Qt 








when Er > Qt, 
when < 



Tr{er) when > Qt- 



(89) 



This case corresponds to 6't ^ 0, or 6*/ ^ 6*0, as can be easily seen from the inspection of Fig. |[ But with this 
geometrical consideration alone, it is hard to find the limit of cosdf in Eq. (|8^). The behavior of the limit of 
cosOf can only be obtained when Fq — s- is taken into account at 9t — > 0, i.e., when the strength of the channel 
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coupling is taken into account, which is hidden in the geometrical consideration because of the use of Eq. When 

-oo when E < Ea 
oo when E > Ea- 

Only two values of Ca are possible in the limit of Fq — s- or \qr \ — > cot A'(2- fa Too correspond to 0q — s- and 
TT, or cos 6*0 1 and —1, respectively. Since cos0f coa 9a, cos Of satisfies the limit ( p8| ) as \qr\ — > cot A^j- 
The energies at which cos 6 f and undergo abrupt changes look different but are equivalent since E = Ea ot 
Ea = corresponds to — Qt owing to the relation ( |75| ) and (ja — > oo. 

3. As IqtI — > oo, the graph of cos^/(er) becomes symmetric with respect to origin and is given by 

cosA?2 

COS 0f[er) — 1^^^=^^^= when \qr \ oo. (91) 
el sin^ A°2 + 1 

The derivation of Eq. ( |9l| ) from Eq. ( ^7|) is not so easy. \qr\ oo arises in two cases, cot A52 — > oo or cosOt = 
{dt — 7r/2). The geometric consideration is a great help when 9t = tt/2. In this case, Napier's rule gives cot 6/ 
= sin^rcot A]'2 [Q- From this formula of 0016*/, cos^Of = cos^ A52/(er sin^ Aj'2 + 1) is obtained by simple 
trigonometric manipulations. The square roots of both sides of it yields Eq. (|9^) except for the sign. In order 
to fix the sign, let us consider the case of = which corresponds to 5r — 7r/2. Since 9t = 5^ = vr/2 means that 
the chord PQ in Fig. ^ is part of the equator, we have Oa = Sa = Tr/2. For this particular spherical triangle, it 
can be easily proved that Of = ■ This fixes the sign. The remaining case of cot A°2 oo corresponds to F^ 
and can not be easily handled by geometric argument as mentioned above. Eq. (74) of Ref. 0] allows us 
to handle this case and yields cosOf = cosAj'j which is identical with Eq. ( ^l| ) in this case. 

The partial delay times at = are 

2n-^(.„^0)^-cos^(^ 

'~dE 

which is easily obtained by substituting A52 for 9f into 2hd6±/dE = Tr{l± cosOf). 

Before doing the data-fitting, let us briefly describe the system used for the calculation and the methods of cal- 
culation. A triatomic van der Waals molecule considered here is of the type of rare gas- • • homonuclear halogen-like 
diatomic molecules [Q. Let us consider the system where the van der Waals molecule in its ground state is excited 
by the light whose energy amounts to the excitation of the diatomic vibronic motion from the u = to 1 state. This 
energy is sufficient to break the van der Waals bond and produces a predissociation spectrum as the energy of light 
is scanned over a certain range of frequency. 

The following interaction potential between A and B2 in A- • • B2 triatomic system 



2ft— (e^ = 0) ^ — sin^ , (92) 



, I4i(i?,r,7) when R < R* 

V{R,r,j)^l „(R-R'Y (93) 

Kdw(r,7) + (Hi-Kdw)e "y^^) when R>R*, 



is the one employed by Halberstadt et. al. to fit the predissociation data for Ne- • • CI2 system and is used here [ p2| . 
In Eq. (p3[), 7 are the Jacobi coordinates that denote the distance between A and the center of mass of B2, the 
bond distance of B2, and the angle between R and r, respectively Vm{R, r, 7) and K-dw are given as 



2 



1=1 



VM(i?,r,7) = DAB}_^i |e-"ABU<AB.-«^BJ _ l _ l ^ (94) 



CM 



(o) 1 2 

-qcm(-R--Rcm' — 1 



(95) 



Kdw(i?,7) = ^6 (96) 
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where RaBi is the distance between A and i^^ B atom; other parameters are adjusted parameters to yield the best 
fit to the experimental values. The values of the parameters used in this paper are given in Table ||. Two Legendre 
terms are retained for Cg{'^) and Cs{j) in Eq. (|9^), e.g., 

Ce{j) = Ceo + Ce2P2icosj). (97) 

R* in Eq. ( p3| ) is chosen as the inflection point of the atom- atom Morse potentials and given by R* — R'^l^ + hi2/ acM- 
The Hamiltonian for the triatomic van der Waals molecules A- • • B2 in the Jacobi coordinates is given in atomic 
units by 

H = -t^tS^ + a + i + V{R, r, 7) + H^^r), (98) 



where 



^^B.(r) = --^|^ + yB.(r), (99) 

denotes the vibrational Hamiltonian of B2; m is the reduced mass of B2; n denotes the reduced mass of A and the 
center of mass of B2; j is the angular momentum operator of B2; 1 is the orbital angular momentum operator of the 
relative motion of A and the center of mass of B2. The values of diatomic molecular parameters of B2 used in this 



paper are given in Table [II . The calculation is limited to zero of the total angular momentum operator J = j + 1, 
as usually done in this field without affecting the predissociation dynamics much. Such a limitation simplifies the 
Hamiltonian as 1 can be replaced by j. 

Let ^'"'•'^(i?, r, 7) be the eigenfunctions of H of Eq. ( p^ ) and let it correspond to the state vibronically excited 
by light which will be predissociating into an atom and a diatomic fragment. It is indexed by the vib-rotational 
quantum numbers {v,j) of its diatomic photo-fragment and abbreviated to i, i.e., i = {v,j). When the wavefunctions 
(R, r, 7) to the dissociation channel i = {v,j) are expanded in terms of n base functions (r, 7) — (r|w')Yj'o(7, 0) 
(«' — 1,2,.. .,n) as 

*-«(i?,r,7)=^$,,(r,7)x,;'.(i?), (100) 
the close-coupling equations for Xi'i{R) ^ire obtained as 

X^'^{R) + J2 V^'^"{R)X^"^{R) = 0, (101) 



1 :2 



2m dR'^ ' 2mi?2 



with 



kl = 2m[E ~ Bj\f + 1) - {v' + i)a>], (102) 

and 

F.'v(i?)= [ d-f sin J [ dr<i>,„{r,j)ViR,r,j)^*,{r,j). (103) 



The close-coupling equation (101) is solved by the De Vogelaere algorithm |34[| and wavefunctions (100) that satisfy 
the incoming wave boundary condition are, consequently, obtained. The S matrix obtained in this process, which 
is identical with (H), is diagonalized and eigenphase shifts (^ are obtained. Two closed channels corresponding to 
{v — l,j — 0) and {v — 1, j = 2) and two open channels corresponding to {v — 0,j = 0) and {v = 0,j = 2) are included 
to mimic the system of one discrete state and two continua to which the theory developed in this work applies. This 
calculation yields the data of eigenphase shifts as functions of energy. Let us call this method of calculation as the 
close-coupling method. 

Note that the theory developed in Ref. and in the present paper relies on the presence of a discrete state 
embedded in continua. Among various theories devised to describe the resonance with explicit consideration of a 
discrete state, Fano's configuration interaction theory is chosen in this work In the normal use as described in 
the above paragraph, the close-coupling method can not be connected with the configuration interaction theory since 
no discrete state is assumed in the close-coupling method. But with a little modification in its use, it can be used 
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to calculate dynamic parameters of the configuration interaction theory. A discrete state with its resonance energy 
£^0 used in the configuration interaction theory can be obtained by solving the close-coupling equations (modified to 
incorporate the shooting method |^ ) with inclusion of closed channels alone. Wavefunctions obtained by solving the 
close-coupling equations with inclusion of open channels alone obviously diagonalize the Hamiltonian in the subspace 

spanned by open channels alone and are the continuum wavefunctions considered in the configuration interaction 

theory. The background scattering matrix 5*" are obtained as a byproduct when the continuum wavefunctions i/'_e^'^ 
forced to satisfy the incoming wave boundary conditions (or outgoing wave boundary condition if a scattering system 
is considered instead of the photo-dissociation). By diagonalizing S*", background eigenphase shifts (5° and ^2 and the 
frame transformation matrix from the asymptotic wavefunctions V'b^''' to background eigen channel wavefunctions 

Ik) . . n . 

■0^ are obtained. With U , background eigenchannel wavefunctions can be obtained from the asymptotic ones as 

^l;f{R,rn) = -^e''"Y.^W\R^rn)^ (104) 
and can be used to calculate partial decay width Tk as 



2-K 



(v4'^|i/|0)|' = 27r|(vi'V(i?,r,7)|0)|', (105) 



where the last equality holds for 5v = ±1 vibronic predissociation since ^(i?, r, 7) is the only term containing odd 
powers of r in H. Since Eq, S^, J^, Ti, and r2 are obtained, the dynamic parameters Ea, Ta, Qa, <Zr, and cot^t are 
calculated using Ea = Eq + AFcot A;2/2, Ta = 2VT\r;/ sin A?2, Eqs. (|), (|^), and (|3|). 

Though the configuration interaction theory directly calculates the dynamic parameters Eq, T, Ea, Ta, Qa, Qt-, and 
cot^t, the close- coupling method can not calculate them directly. Dynamical parameters directly obtainable from 
the close-coupling method are the S matrix, its eigenphase shifts 5+ and 5-, and partial delay times 2hdS+/dE and 
2hdS-/dE as functions of energy. If the assumptions used in configuration interaction theory are exact, eigenphase 
shifts and partial delay times should satisfy Eq. (^) and Eq. (^7|), respectively. The assumptions that S^, S2, Ti, and 
r2 are constant of energy on which the configuration interaction theory relies are expected not to cause much trouble 
in the actual situation as far as the energy range is not wide enough. If eigenphase shifts and partial delay times 
follow Eq. (H) and Eq. (p7|), values of dynamic parameters Ea, Ta, Qa, Qt, and cot^t can be obtained by varying 
them so that eigenphase shifts and partial delay times calculated by the close-coupling method fit the formulas best. 
Fitting is done with the Levenberg-Marquardt method recommended for the nonlinear models in Ref. p5| . For the 
reason mentioned above, only partial delay times will be fitted. 

Data fittings are done in two steps. At first, Eo and F are obtained by fitting the eigentime delay sum. Then q-r 
and cot6't are obtained by fitting partial delay times 2hdS+/dE {=Tr + Ta) and 2hd6-/dE (—Tr — Ta) with (p7|). Either 
2hd6+/dE or 2hd6-/dE can be used to obtain the values of qr and cot 6*4. Values of g,- and cot 6't obtained from either 
of them should be identical if the assumption of the configuration interaction method, namely that Fi, F2, S^, and S2 
are constant of energy, is exact. The differences between qr's and cot 6't 's for 2hdd+/dE and 2hd6-/dE may serve as 
a criterion of the exactness of the configuration interaction theory. 

Numerical study shows that fitting the eigentime delay sum calculated by the close-coupling method to the formula 
( p^ ) is done reliably. Compared with that, fitting the partial delay times to the formula ( ^7|) can not be easily done. 
Eigentime delays calculated by the close-coupling method show abnormal behaviors like negative values of eigentime 
delays at the energy region not far from a resonance while theory contends that eigentime delays are positive in the 
neighborhood of a resonance, meaning that assumptions made for the formula ( p7| ) are prone to break down. It may 
be argued that parameters obtained by data-fitting of eigentime delays calculated by the close-coupling method will 
deviate more from those obtained by the configuration interaction method as the range of energy taken for fitting 
is wider. But we can not always narrow the range of energy for this reason. Notice that qr is the energy where 
meets r_. As seen in Eq. (|89|), if qr is large, two curves are almost identical with except for the neighborhood of 
Qr- Therefore for the good fitting, the energy range taken for fitting should include qr- It means that if qr is large, 
fitting is likely to be bad. This assertion is checked by doing data-fittings for three different ranges of energy, namely, 
[E^ - F, + F], [Eo - 2T, Eo + 2T], and [£^0 - 3F, Eo + 3F]. 

Table IV shows the results calculated with the parameters given in Table || and |l|. Table is obtained with 



the same parameters but with — 5.044 a.u. The former table corresponds to the case of large qr while the latter 
table to the case of small qr- The tables also show values of a fudge factor A which can be used as a criteria for the 
goodness-of-fit. The fitting is done at first by steepest decent method with the initial value, say 0.001, of A. New 
value of A is suggested for the next iteration. If the suggested value of A becomes sufficiently small, inverse-Hessian 
method is used for fitting. At the final call, A is set to zero. This method assumes that values of A should approach 
zero if everything goes O.K. This is true with the data-fitting of the eigentime delay sum. For the partial delay times. 
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values of A do not go to zero. Though values of A do not go to zero, comparison of curves obtained by two method 
showed that the fitting is good enough to be acceptable if values of A are not too large. In Table the values of 
qr lies beyond the first interval [Eo — F, _Eo + T] and lies in the second interval. According to the reasons mentioned 
above, it is hard to achieve reliable data-fitting in this case. The closest result to the theoretical values is obtained 
for the first interval, the narrowest one, where the fudge factor is worst indicating wider interval should be used. The 
calculation shows that wider interval yields worse result indicating that the assumptions may no longer be true. On 
the other hand, in Table V, the value of qr is small and the data-fitting may be done rather reliably and is confirmed 
by the calculation. This situation contrasts greatly to the fitting of partial photo-dissociation cross-sections to the 
Fano-Beutler-like line profile formulas [E9[ 



<Jj=^° \ , = ^1 9 + < . ? ' (106) 



for the same predissociating system of van der Waals molecules, where the fitting is excellent as shown in Tabic VI 



VII. SUMMARY AND DISCUSSION 

In the previous work [Q, eigenphase shifts for the S matrix and Smith's lifetime matrix Q near a resonance were 
expressed as functionals of the Beutler-Fano formulas using appropriate dimensionless energy units and line profile 
indices. Parameters responsible for the avoided crossing of eigenphase shifts and eigentime delays and the change 
in frame transformation in the eigentime delays were identified. The geometrical realization of those dynamical 
parameters is tried in this work, which allows us to give a geometrical derivation of the Beutler-Fano formulas 
appearing in eigenphase shifts and time delays. 

The geometrical realization is based on the real three-dimensional space spanned by the Pauli matrices fJx, o'y, cr^ 
as basic vectors, where vectors are orthogonal in the sense that 

Tiia.aj) ^2S,j. (107) 

Such a kind of space is called a Liouville space . A 2 x 2 traceless Hermitian matrix, which is generally expressed as 
(T • r with r real, is a vector in this Liouville space. The magnitude of a vector in the three-dimensional Liouville space 
corresponds to the degree of anisotropy in coupling of eigenchannels with the cause that brings about the dynamics 
of the dynamic operator corresponding to the vector. The four-dimensional Liouville space including the unit matrix 
1 as another basic vector is also used for the 2x2 Hermitian matrix whose trace is not zero. 

Resonant scattering can be separated from the background scattering in the S matrix for the multichannel system 

around an isolated resonance as 5 = 5°(7rfc + e^'^^^''Tra) where tTq is the projection matrix to ipE^ which is the only 
type of continua interacting with the discrete state and tt^ = 1 — tt^ [pO| . When the number of open channels are 
limited to two, 5°, ni, + e~^^^''TTa, and S can be expressed using Pauli's spin matrices as 

^ ; 
g-i(5El+i5j,<T-na) 

Phase shift matrices A , A^, and A which are Hermitian may be defined for and S as 

A° = i(<5°l + AV-^), 

Ar = ^((^rl + SrO- ■ fit), 
A= ^{5^1 + 6a(T ■ Ua) 

and are vectors in the four-dimensional Liouville space. 

According to the Campbell-Baker-Hausdorff formula, the phase shift matrix A is expressed as an infinite sum of 
multiple commutators of A" and A^, which is difficult to use |2^]. The geometric way of representing the combining 
rule of a" and A,, into A provides an alternative to that. We first note that the isotropic part of A is simply 
obtained from those of A° and A,, as a simple addition -I- Sr) and factored out in iS = iS'^(7rf, -|- e~^^^^''TTa)- Then 
the remaining anisotropic part of A is obtained from those of A*^ and A^ as 



s = 
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The above identity can be expressed using the rotation matrices in the Liouville space as 

This relation leads to the construction of the spherical triangle whose vertices are the endpoints of vectors correspond- 
ing to the anisotropic parts of the phase shift matrices but whose lengths are limited to unity. The lengths of the 
vectors are utilized as the vertex angles of the spherical triangle. This spherical triangle shows the rule of combining 
the channel-channel couplings in the background scattering with the resonant interaction to give the avoided crossing 
interactions in the curves of eigenphase shifts as functions of energy. 

The time delay matrix Q basically derives from the energy derivative of the phase shift matrix A of the S matrix. 
The phase shift matrix A is a vector in the four-dimensional Liouville space and is ^{Ssl + 6acr ■ fia) as stated above. 
The energy derivative of 5^ and 5a without changing the direction of the vector fia yields the "partial delay time 
matrix" given by h{d5x-^/ dE + cr ■ ha)dda/dE. Discussion on the time delay matrix due to a change in the direction 
of the vector fia or in frame transformation is greatly facilitated by the use of the formula ^hidO a / dE)[S'^ cr ■ yS — 
cr ■ y). The formula shows that the time delay due to the change in frame transformation is the interference of two 
terms. The first term inside the parenthesis comes from the energy derivative of the frame transformation from the 
background eigenchannels to the S matrix eigenchannels and the second term from the energy derivative of the frame 
transformation from the S matrix eigenchannels to the background eigenchannels. The first term corresponds to a 
vector rotated from the y vector by the rotation Rn^ {~25a)- The net time delay resulted from the interference is thus 
calculated by the vector addition in the three-dimensional Liouville space. 

Going back to the spherical triangle, the laws of sines and cosines and other laws derivable from such laws such as 
the laws of cotangents holding for the spherical triangle can be translated into dynamic laws by converting 5r and 9a 
into energies according to cot^^ = —^r and cot^a = — e^- Two cotangent laws 

sin 1^12 cot 5a = — sin 9a cot 9t + cos 9a cos , 
cot 0/ sin(?t = — cos(7r — 5r) cos 9t + sin(7r — 5r) cot AI2, 

can be shown to correspond to two Beutler-Fano formulas for cot 5 a and cot6'/ 

— Qa 



cot 5a — ~ cot 



cot9f = — cot I 



Er - Qr 



with such conversion. Other laws also yields interesting relations among dynamical parameters. 

The presence of the dual triangle of the spherical triangle indicates that we can make a one-to-one correspondence 
between edge angles with vertex angles so that if there is one valid relation we can make another valid relation 
by interchanging angles with its one-to-one corresponding angles. In other words, for each edge angle, we have a 
conjugate vertex angle and vice versa. This conjugation relation can be extended to and by making use of their 
relation with 5r and 9a, respectively. The full conjugation relations among geometrical and dynamical parameters are 
listed Table ||. The duality of the spherical triangle thus explains the symmetry found in the dynamic relations and 
provides us with a systematic approach and complete symmetric relations. Besides this use of trigonometric laws of 
the spherical triangle, the geometric construction in the Liouville space facilitates other useful consideration. 

Note that the geometrical laws holding for the geometrical objects in the real three-dimensional Liouville space deal 
only with the intrinsic nature of the dynamic couplings independent of the characteristics of the individual system. 
It derives from that the reduced energies hide the specific characteristics of the dynamic couplings of the individual 
system such as the strengths of the dynamic couplings between the discrete state and continua, the indirect couplings 
between continua via discrete states, the resonance positions, the avoided crossing point energy. Intrinsic nature of the 
dynamic couplings is concerned with the relations among eigenchannels for various dynamic operators, the anisotropy 
in the channels coupling in the 5"° and S scattering, the anisotropy in the channel coupling with the discrete state. 
This shows both the beauty and the limitations of the geometrical construction in the Liouville space. In the actual 
application, we have to be careful when considering the case close to the limits in coupling strength, where abnormal 
behaviors take place in actual dynamic quantities, but where no abnormality shows up in the Liouville space. 

The present theory is developed for the system of one discrete state and two continua. It will be desirable to extend 



the theory to more than two open channels and to overlapping resonances for which the results of Rcfs. |17 19[| will 
be a great help. It might be also valuable to apply the present theory to MQDT. In connection with the latter, it 
might be interesting to apply the present theory to extend the Lu-Fano plot to multi-open channel case. 
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APPENDIX A: SEPARATION OF THE RESONANT CONTRIBUTION FROM THE BACKGROUND 

ONE IN S MATRIX 



The S matrix for the multichannel system in the neighborhood of an isolated resonance given in Eq. (^ can also 
be written as 



VieV, 



kE 



E-Eo- ir/2 



(Al) 



where Vie is defined in ter ms o f the asymptotic state V's*'"''* (''Pe'^'" l^l^*)- Using the relation (E — Eq + iT/2)/ (E 



Eo - ir/2) 



Eq. (Al) becomes 



-2iSr 



1) y^ieV, 



kE 



(A2) 



If we consider Fano's ^/j^ ' state defined as 



T.\^EnVjE 



(A3) 



and the projection operator |V'b'')(V'e'*I into that state, and define the projection matrix Ha whose (Z,j) element is 
given by 



(riah = {^E^'^\^^E^){^'e^\^e'") = =^V,eV;, 



27r. 



Eq. (A2) becomes in matrix form as 



-2iSr 



l)5"n„ = 5" (Hb 



(A4) 



(A5) 



where H;, is equal to 1 — Ila and is the projection matrix into th e sp ace orthogonal to the space of ?/'^'' and may be 
regarded as a projection matrix into the background. Using Eq. (A5), the time delay matrix Q becomes 



Q = 2ft^n,. 

dE 



(A6) 



Notice that Fano's tA^"* and other continua orthogonal to it are the eigenstates of the time delay matrix Q for the 
multichannel system in the neighborhood of an isolated resonance. 

Now let us confine the number of op en c hanne ls of the system to 2. Then the background scattering matrix S''^ and 
projection matrices and H;, in Eq. (A5) and ( |A6| ) can be expressible in terms of Pauli matrices. When we convert 
them into ones in terms of Pauli matrices, it may be convenient to choose the background eigenstates as a basis. Let 
us denote the background scattering matrix and projection operators in the basis of background eigenstates as , 
TTa, and TTf,. They can be obtained from those in the basis of asymptotic channel wavefunctions as U°S'^U''\ U'^HaU'^, 
U^IlbU'^, respectively. Eqs. (|A5|) and (|A6|) are now in the basis of background eigenchannels as follows 



S = S° (tT, + e-'^'^TTa) 

Q = 2h^na. 
dE 



(AT) 
(A8) 
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Now in the basis of the background eigenchannels the background scattering matrix 5*^ takes the form 



5' 



Comparing Eq. (p9[) and Eq. (A8), we obtain 



e-^'^i 



= e 



-i(<5^1+Ai'2'T-«) 



(A9) 



from which TTf, is obtained as 



TTa = - (1 + (T • fit) 



1 , 

nt ^ -(1 - (T ■ Tit) . 



Inserting Eqs. ( AlC ) and ( All ) into Eq. (A7), we obtain 



(AlO) 



(All) 



(A12) 



Equating Eqs. (|5§) and ( |A12|) , Eq. (|6g) is obtained 



APPENDIX B: MULTIPLICATION OF TWO ROTATIONS AND SPHERICAL TRIANGLE 



Multiphcation of two rotations yields another rotation. The question is how the three rotations are related geo- 
metrically in real three-dimensional space. If three rotations are represented by three rotation axes fiA-, fiB, nc and 
corresponding three rotation angles 2A, 2B, 2tt — 2C, we found that three rotation axis form a spherical triangle 
whose vertex angles are given by A, B, and C (the corresponding opposite edge angles will be denoted as a, b, and c). 
When this theorem is expressed in terms of unitary matrices in the complex two dimensional space, it takes the form 



(Bl) 



The proof of this theorem can easily be done if we use the relation e-*<"^ " = cos a—icr-h sin a and various sine and cosine 
laws of the spherical triangle and will be omitted. The similarity transformation g-^i^-C)<T-nc j^^i('^-C)eT-nc jg equal 
to e-^BeT-nB^-iAeT nAj^giAfT nAg.tBrT nB owing to the identity (Bl). Such an equality of two similarity transformation 
can be expressed in terms of rotation matrices as 



RnA{2A)Rn^{2B) = Rnc{~2C). 



(B2) 



Sometimes, it might be convenient if rotations in the identity (Bl) are expressed in terms of basic axes instead of 
arbitrary ones. Without loosing generality, let us take fiA as the z axis and let hb lie on the zx plane. Then the 
spherical polar coordinates of three vectors are ua — (1,0,0), ub = (l,c,0), and nc — (1,^,^)- Using Rng{2B) = 
Ry{c)Rz{2B)Ry{—c) and similar formula for _Rft^(— 2C), Eq. (|B2| ) becomes 

RrCA{'^A)Ry{c)Ri{2B)Ry{^c) = Ri{A)Ry{h)R,{-2C)Ry{-h)Ri{-A). (B3) 



The validity of Eq. (B3) can be proved by inserting the explicit formula for the rotation matrices, for example, such 
as 



RziO) 



cos 6 — sin 9 
sin6' cos6' 
1 



(B4) 



and by using sine and cosine laws of spherical triangle, the derivation of which is straightforward but requires patience. 
The equality of two relations (B2) and (B3) can be shown using the finite rotation formula for the vector r' which is 
obtained by rotating r about the unit vector fi through an angle $ as follows 



r' = n{n • r) + [r — n{n ■ r)] cos<I> + (r x n) sin$. 



(B5) 
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TABLE I. Conjugate variables. 



Variable Conjugate variable 

5r — Oa 

Of 5a 

A% n-9t 

€BF,r CBF.a 

qr -go 



TABLE IL Parameters for the model intermolecular system A- • • B2. 

Parameter Value 

Reduced mass between A and B2 

m 6756.8 a.u. 

Morse potential parameters 

Dab 0.0034 eV 

q;ab 1.0 a.u.~^ 

Ri°^ 6.82 a.u. 

DcM 0.00195 eV 

acM 1-0 a.u.~^ 

R^°2i 6.65 a.u. 

Van der Waals potential parameters 

Ceo 0.75 eV(a.u.)-^ 

C62 0.119 eV(a.u.)-« 

Cso 1.58 eV(a.u.)-* 

C82 0.8 eV(a.u.)-* 



TABLE III. Diatomic molecular parameters. 

Parameter Value 

vibrational frequency We 0.0162 eV 

rotational constant B 0.01758 meV 

equilibrium bond length Ve 3.044 a.u. 

reduced mass fx 32576.6 a.u. 
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TABLE IV. Comparison of parameters obtained from the theory in this paper and from the data-fitting of partial delay 
times calculated by the close-coupling method. The channels included in the close-coupling equation are described in the text. 
Values of parameters in Tables O and III are used in the calculation. The three energy intervals used for the data-fitting in 
the close-coupling method correspond to [—2,2], [—4,4], and [—6,6], respectively, in scale. In this system, cot^ dt 3> 1 and 
9t ~ — TT. This means that q^- — cot A12, corresponding to the limiting case considered in Eq. (|89|). Note that lies in the 
second interval [—4, —2] beyond the first interval. The avoided crossing point energy Ea is located at —4.4 x 10~^ eV from the 
resonance energy Eq and is very close to q-r when measured in scale since it is equal to qr cos^ 9t . sin 9t ~ means that r2 ~ 
0, physically corresponding to the system where the decay process into one open channel is much faster than that into another 
open channel. 





Theory 




Close-coupling method 






in this paper 


[Eo - r, £0 + r] 


[Eo - 2r, Eo + 2r] 


[Eo - 3r, Eo + 3r] 


£;o(eV) 


0.01396397 


0.01396567 


0.01396567 


0.01396567 


r(eV) 

It 


2.38x10"^ 
-3.69 


2.35x10"^ 
-3.56 


2.34x10-^ 
-3.25 


2.320x10-^ 
-3.45 


cot 2 Ot 


5274 


5000 


1.0x10^ 


2.7x10^^ 


A 

qr 


-3.69 


10-1 
-3.56 


lo-'^ 

-3.40 


10-* 
-3.51 


COt^ dt 


5274 


5000 


5.9x10* 


-1.4x10* 


A 




10-1 


10-^ 


10-^ 



TABLE V. Comparison of parameters obtained from the theory in this paper and from the data-fitting of the partial delay 
times calculated by the close coupling method. The channels included in the close-coupling equation and values of parameters 
used in the calculation are the same as those in Table ^ except that re = 5.044a.u. Avoided crossing point energy is located 
at 8.1 X 10-*eV from the resonance energy and lies very close to the resonance energy, meaning that A?2 ~ 7r/2. Then qr ~ 0. 





Theory 




Close-coupling method 






in this paper 


[Eo -r,Eo + r] 


[Eo - 2r, Eo + 2r] 


[Eo - 3r, Eo + 3r] 


£;o(eV) 


0.01553487 


0.01553486 


0.01553486 


0.01553486 


r(eV) 

q-r 


3.35x10-^ 
0.115 


3.29x10"^ 
0.164 


3.27x10-*^ 
0.199 


3.24x10-*^ 
0.254 


cot 2 dt 


0.73 


0.71 


0.75 


0.78 


A 

qr 


0.115 


10-* 
0.109 


lo-'^ 

0.132 


10-^ 
0.133 


COt^ dt 


0.73 


0.71 


0.74 


0.77 


A 




10"=* 


10-^ 


10-^ 



TABLE VI. Comparison between line profile indices of partial photo-dissociation cross sections obtained by the configuration 
theory and by the close-coupling method. Adopted from Ref. [29]. 



j 


Confif 
(T° [arb. unit] 


juration interaction method 


mq,)? 


Close-coupling method 
(7° [arb. unit] ^(qj) 


mq,)? 





0.221 


-356 


24280 


0.219 


-358 


24920 


2 


0.279 


-233 


169 


0.281 


-233 


-48 


4 


0.305 


-171 


23810 


0.310 


-173 


23040 


6 


0.282 


-275 


10930 


0.281 


-275 


11390 


8 


0.154 


-248 


34190 


0.150 


-248 


36880 


10 


0.035 


143 


76130 


0.034 


125 


81320 
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FIG. 1. Interference between two time delay processes due to the change in frame transformation. 



FIG. 2. Four coordinate systems xyz, x'y'z', x"y"z", and x"'y"'z"' pertaining to the eigenchannels of iS", S or a Va, crVf, 
and a Vt- 



FIG. 3. A diagram showing the spherical polar coordinates of fit in the xyz and x'y'z' coordinate systems. The projection of 
ht on the x'y' plane coincides with hf except for the length. Note that ha, fit, and fif lie on the same plane. This diagram shows 
that the spherical polar coordinates of ht in two coordinate systems are given by {l,dt, —^12) and (1,^/, —Sa), respectively. 



FIG. 4. The spherical triangle whose vertices are the endpoints of z, ha, and ht- 

FIG. 5. The dual spherical triangle corresponding to the spherical triangle of Fig. ^, where the edge angles of the latter 
becomes the vertex angles of the former and vice versa. 

FIG. 6. A diagram showing the geometrical traversal of the vertex P as energy increases from —00 to 00, the point P traverses 
from the point A to the opposite point —A along the great circle. 



FIG. 7. Partial delay times 2hdS±/d€r vs. er are plotted for three different profile indices ~ 0.6, 1, and 5 with A12 = tv/3. 
Graphs for 2hdS+/der are located more shifted to the left hand side than the corresponding ones for 2hdS- / der- The values 
0.6 of is close to its possible minimum value 0.577 (=cot Aja). 
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